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ABSTRACT 

> ■ 

• i-H , 

We present the evolution of the color-magnitude distribution of galaxy clusters 
from z = 0.45 to z = 0.9 using a homogeneously selected sample of ~ 1000 clus- 
ters drawn from the Red-Sequence Cluster Survey (RCS). The red fraction of galaxies 
decreases as a function of increasing redshift for all cluster-centric radii, consistent 
with the Butcher-Oemler effect, and suggesting that the cluster blue population may 
be identified with newly infalling galaxies. We also find that the red fraction at the 
core has a shallower evolution compared with that at the cluster outskirts. Detailed 
examination of the color distribution of blue galaxies suggests that they have colors 
consistent with normal spirals and may redden slightly with time. Galaxies of starburst 
spectral type contribute less than 5% of the increase in the blue population at high red- 
shift, implying that the observed Butcher-Oemler effect is not caused by a unobscured 
starbursts, but is more consistent with a normal coeval field population. 
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Subject headings: methods: statistical - galaxies: elliptical and lenticular - galaxies: 
evolution - galaxies: clusters: general 



Introduction 



The comparison of local and high redshift clusters of galaxies is a powerful tool for constrain- 
ing theories of structure formation and cosmological world models. The discovery that galaxy clus- 
ters at z ~ 0. 5 have a larger fraction of blue g alaxies compared with similar clusters found in the lo- 
cal universe (But cher & Oemlenll978Ul984|) provided the first direct evidence of galaxy evolution 
in dense environ ments. Since these early st u dies, the Butcher-Oemler effect has been confirmed 
photometrically (Rakos & Schombert 1995 ^ Margoniner & de Carvalho 200ol: Margoniner et al 



2001 



Kodama & Bowerl 



1999 



Ellingson et al 



2001 



Goto et al.l 



2001 



2003b and spectroscopically ( Dressier & Gunnl 19821: 



Laverv & Henrvll988:Wbricant. McClintock & Baud 199 lclbressler & Gunnlll992l : IPoggianti et al. 



Pogg ianti et al. This observed evolution has been interpreted 



variously as evidenc e for a rapid change in galaxy populations driven by ram-p ressure stripping 



(Gunn&Gott 



1988 



Moore et al 



19721). ga laxv-merging or other physical processes in clusters (|Lavery & Henry 



1996), as well as the effect of infalling galaxies from the field or near- field 



regions ^auffmannl 1995 : Abraham et al. 1996 : Ellingson et al. 2001 ) 



However, the Butcher-Oemler phenomenon has been challenged from the start (e.g., iKoo 



198 If) and more recently by a host of new analyses, dra wn primarily from X-ray selected sam- 
ples of clusters and near-infrared observa tion of galaxies (|Andreonlll998l : lAndreon & Ettorill 1999c 
Smail et al.lll998uDe Propris et alj|2003j) . In these, the evolutionary trend with redshift was found 
to be either weaker, or non-existent. The non-uniformity of the original But cher-Oemler clusters, 



espe cially in physically relate d properties like galaxy surface number density (lNewberry. Kirshner & Boroson 



19881) and X-ray luminosity (lAndreon & Ettorilll999|) may present selection biases that lead to a 
strong redshift evolutionary trend. Indeed, the original Butcher-Oemler clusters were selected 
from photographic plates and may include lower mass and unvirialized merging systems. The 
strong ^-correction pushes the sample to bluer rest bands. These opticall y-selected galaxy sam- 
ples may also prefere ntially select systems with strong star-formation ("e.g. jGunn. Hoessel & Oke 
19861 : lAndreonll2005l) . 



Earlier studies also suffered from low number statistics. While the analyses were usually done 



'based on observations carried out at the Canada-France-Hawaii Telescope (CFHT), operated by the National 
Research Council of Canada, the Centre National de la Recherche Scientifique of France and the University of Hawaii. 

2 based on observations taken at the Cerro Tololo Inter- American Observatory 
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with fewer than 50 clusters over a large lookback time (e.g., from z = to z = 0.5), general claims 
were often made about the cluster population as a whole. Part of the reason is that investigators 
necessarily analysed whatever systems that were available. Uniformly selected cluster samples that 
span a large redshift range and have constant physical attributes (e.g., similar m asses) are sparse 
and th e variation of blue fraction within a fixed epoch is often considerable (e.g jMargoniner et al. 



200 II) . It is unclear to what extent this variation is due to actual cosmic variance or the perculiarities 
of the sample selection itself, since the number of clusters involved is at most a few tens. 

With the advent of wide area ph otometric and spectroscopic redshift surveys like the Sloan 



Digital Sky S urvey (I York et all 12000 . SDSS) and the Two-Degree Field Galaxy Redshift Survey 



at low (z ~ 0.1) to intermediate (z ~ 0.3) redshift 


'e.g.. 


Annis et al. 


1999; 


Goto et al. 


2002 




Kim et al.ll2002 


: Merchan & ZandivareJ2002 


;l 


Eke et al.l2004 


Yang et al. 


2005 


;l 


Merchan & Zandivarez 


2005; 


Weinmann et al. 


2006; 


Miller et al. 


200f 


i;lBerlind et al. 


2006 


). However, in order to under- 



stand what processes are dominant in explaining the Butcher-Oemler effect, a large redshift range 
is desirable to provide maximal leverage to evolutionary studies as various theoretical models may 
have different redshift dependences. 

In this work, we revisit the Butcher-Oemler phenomenon using the currently largest sam- 
ple of uniformly selected clusters of galaxies available, out to z ~ 0.9. The ga laxy clusters are 
drawn from the complete and well-calibrated Red-Sequence Cluster Survey (RCS: lGladders & Yee 
2005J). The RCS method searches for galaxy overdensities along a color-magnitude sequence cor- 
responding to cluster early-type galaxies, the cluster "red sequence". The apparent color of this 
sequence reddens with redshift, allowing us to separate overlapping structures and effectively re- 
move most of the contamination that has plagued monochromatic cluster surveys. All massive 
clusters discovered by any methodology to date show a detectable red sequence, so that, while it 
is slightly less effective for very blue systems, this methodology is not expected to introduce a 
strong bias to the study of the cluster populations. To the extent that galaxies on the early-typ e 
red sequence have the largest fraction of stellar mass at a fixed luminosities (e.g. jBell et a l. 2004), 
our cluster sample can be considered stellar mass dependent. Our sample may exhibit selection 
differences when compared, for example, with d ark mass dependent weak lensing detected cluster 



samples such as the Deep Lens Survey (DLS; Wittman et al 



20021) . intracluster medium (ICM) 



density dependent X-ray selected samples (e.g..: iRosati et al] |1998|) . or star-forming galaxy bi- 
ased blue overdensities for Lyman break selected samples. However, galaxy density maps derived 
from red seq uence galaxies have been shown t o correlate with weak lensing shear maps with high 
significance (jwilson. Kaiser & Luppino 2001 ). 



In this paper we study the evolutionary trend of the red fraction for cluster galaxies with red- 
shift within a range of radii, scaled by r/r 2 oo - the radius at which the galaxy overdensity is 200 
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times the mean density. Scaling the radii to r/ r 2 oo allows a statistical combination of many systems 
and minimizes many problems arising from t he foreground/background subtraction in photometric 
surveys. Our approach is clu ster-centric (e.g jMemick & Sargenq|l977l : IWhitmore. Gilmore & Jones 
19931 : LEUingson et al.ll2001|) . where we use a characteristic radius r 2 oo to normalize the clusters 
across the richness range of the clusters. This approach is complementary to many more recent 



studies of the C olor (Star Formation RateVDensit v relation (IHogg et al 



Yee et al 



2005 



2005 ) which use local galaxy densitie s (IDresslenll980l : ISmith et al 



Cucciati et al.ll2006l : ICooper et al.ll2007al) . 



2003 



Balogh et al. 2004; 



2005 



Postman et al. 



The outline of this paper is as follows. In the next section, we describe the sample used in our 
analysis. The procedure for the statistical background correction of composite clusters is described 
in section[3] In sectionlU we discuss our results and give caveats that need to be taken into account 
for robust interpretation. In section [51 we summarize our conclusions. Throughout this paper, 
h = ifo/100 km s" 1 Mpc -1 and we adopt the concordance Vt m = 0.3, = 0.7 cosmology. 



2. Sample 



2.1. The RCS Survey 



The Red-Sequence Cluster Survey (RCS; iGladders & Yee! 120051) is a ~ 100 deg 2 imaging 
survey in the z' and R c bandpasses conducted with the CFHT 3.6 m and CTIO 4 m with the 
primary goal of searching for clusters of galaxies up to z fa 1.4. There are 22 RCS patches 
(each of ~ 2.5 deg 2 ) distributed over the northern and southern sky. Some of the patch es were 



ch osen to overlap with o ther widely studied fields, s uch as XM M-LSS (IPierre et al 



20041) . CDF- 



S (IGiacconi et all 120011) . CNOC2 (lYee et all 1200 11) . SWIRE (lLonsdale et all 120031) and PDCS 
IPostmanetal.il 1996k 



2.2. Photometry, Calibration, Extinction 



The primary RCS photometric reduction of object detection, source deblending, star-galaxy 
separation and aperture photometry w as performed using the Picture Processing Package (PPP; 
Yedll99ll : lYee. Ellingson & Carlberdl 19961). Th e detailed implementation of the RCS photometric 
pipeline is described in IGladders & Yed (I2005|) . The average 5cr limiting magnitudes for point 
sources are z' = 23.9 and R c = 25.0 (AB system). Analysis of log(N)-log(S) galaxy counts 
suggests that the photometry for extended sources is complete to ~ 0.8 mag brighter, although 
there are variations from pointing to pointing. There are over 3 million point sources and 11 
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million gala xies in the RCS d atabase. The relative photometry between patches is better than 0.1 
magnitudes (Hsieh et al.ll2005|) . Each d etected extended source is further c orrected for Galactic 
extinction, interpolating the dust map of lSchlegel. Finkbeiner & Davisl (I1998|) and assuming R Rc = 
2.634 and R ,i = 1.479. For the present analysis, we use 72 deg 2 of the RCS survey with the best 
photometry dGladders et al]|2007l) . Throughout this paper, our R c and z' photometry is presented 
on the AB system. 



2.3. The Cluster Sample 



About 8000 groups and clusters of galaxies were detected with high significance by search- 
ing for source overdensities in color slices along th e red-sequence of early -type galaxies. De- 
tails of the cluster detection method are outlined in IGladders & Yed (12000). The release of the 
first catalog, as well a s improvements over the cluster detection methodology, are described in 
Gladders & Yed (I2005|) . A useful by-product of the red-sequence method is an estimated photo- 
metric redshift for each detected cluster. RCS uses only two bandpasses strategically chosen to 
straddle the 4000A break of the target redshift covera ge, and photometric redshifts of the detected 
clusters are quite well constrained, to Az/z pa 0.1 (IGilbank et al.ll2007|) . The dense collection 
of red galaxies along the red-sequen ce effectively diminishes the Poisso n noise from photometric 
uncertainties of individual sources (ILopez-Cruz. Barkhouse & Yedl2004) . 



In addition to an estimated photometric redshift (hereafter z ph ), an o ptical richness parametrize d 
by the amplitude of the galaxy-cluster cross- correlation function B gc (|Longair & Seldnerill979 : 



Yee & Greenl Il987l : lYee & Lopez-Cruzl 119991) is measured. B gc is found to correlate well with 



cluster properties such as X-ray temperature and galaxy velocity dispersion at intermediate red- 
shift, and it is relatively well calibrated for X-ray selected clusters at z ~ 0.35 (Yee & Ellingson 
2003 ). Note that here we use red sequence galaxies for the measurement of B gc (|Gladders & Yee 
2005 1), which give s a tighter relation with respect to X-ray properties, especially at higher redshift 
(|Hicks et al.ll2007|) . ROSAT observations have also found that most g roups and lower mass sys- 
tems with extended X -ray halos are dominated by early-type galaxies (|Osmond & Ponmanll2004t 
Mulchaey et al.ll2003|) . B gc computed from the red sequence will down-weight unvirialized poor 
systems without a hot ICM. 

For the present analysis, we use ~ 1000 clusters with B gc > 300 (h^ Mpc) 177 and 0.45 < 
2ph < 0.90. Mass calibration using low to intermediate redshift clusters suggests these sys- 
tems have richn esses corresponding to Abell Richness Class (A RC) > 0, o v > 450 km s™ 1 and 
Tx > 3 keV (|Yee & Lopez-Cruzl 1 19991 : lYee & Ellingsonll2003l) . The sample is divided into five 
redshift slices with 171 to 263 systems in each slice (Table [l]). While the mass of each cluster 
is uncertain due to the rather large uncertainty in the measured B gc and the intrinsic scatter in 
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mass-£> gc relation, the overall distribution of cluster counts as a function of redshift is robust. Fig- 
ure Q] shows the redshift histogram of clusters used in our analysis. Analysis of the RCS cluster 
counts to constrain co smological parameters yields results consistent with the concordance values 
( Gladders et al. 2007|) . indicating that large systematic errors in our mass estimates are unlikely. 
From abundance matching, clusters with B gc > 300 (h^Q Mpc) 177 used in this paper are rare, with 
comoving densities < 10~ 5 /i 3 Mpc~ 3 , consistent with dark matter halos with mass > 10 14 /i 1 M . 



3. Methodology 
3.1. Foreground and Background Subtraction 



To obtain the composite color-magnitude distribution of our cluster galaxies, we have to sub- 
tract the contribution from foreground and background galaxies statistically. The two-dimensional 
observed background counts for non-stellar sources in R c — z' color and z' magnitude are con- 
structed for each of the 22 RCS patches. Depending on the relative number of clusters from each 
patch used for the composite cluster analysis at each redshift slice, a weighted background based 
on the effective area contributed by the clusters is created (Figure [2]). This semi-global background 
approach allows us to model the background more accurately using a large area, but still mimics 
some inherent patch t o patch variations s uch as observing depth and seeing. This approach was 



z 



Hansen et al 



(2005) but performed in two dimens ions (in color and magni 



similarly employ ed by 

tude s pace) as in LLoh & Strauss! (120061) . O ur approach is in contrast with Valotto, Moo re & Lambas 



(2001) who advocate a local background subtraction methodology. iGoto et al. d2003h find little 
difference with the two approaches if a large enough area is used. The global approach has the 
advantage that it uses all available data to constrain a smooth background with high precision. For 
a given patch of survey area, the region occupied by the clusters, after aggregating in redshift, takes 
up 2 — 5% of the entire patch area. We rescaled the background to match this total cluster area 
before the background subtraction. 

We are interested in the average properties of clusters at a given redshift. If we assume 
that clusters are self-similar systems with a universal galaxy density profile, then clusters of dif- 
ferent physical sizes can be scaled to some normalized cluster radius. However, without X- 
ray information or spectroscopy, it is difficult to measure the virial radius of each of the clus- 
ters. We use the measured B gc for each cluster to estimate a virial radius r 2 oo using the relation 
logr2oo = 0-48 log B zc — 1.10. This relation is empirically calibrated from intermediate redshift 



clusters (|Yee & Ellingsonll2003l : lBlindertl l2006). We stack all clusters in a given redshift slice, 
but scale each cluster by its r 2 oo for an effective normalizing scale. We emphasize here that B gc 
is measured at the core, within 0.25 h^ 1 Mpc, and that the self-similar assumption is implicit in 
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Yee & Ellingsonl(|2003|) . We note here that the spherical halo assumption inherent in the B gc mea- 



surements is not perfect, as in addition to spherical distribution of red galaxies, the RCS method 
may systematically include physically associated filamentary galaxy structure along the line-of- 
sight. However, because of the overdensities implied at B gc > 300 (h^ Mpc) 177 , most structures 
will be dynamically collapsed, if not completely virialized. We are exploring secondary measure- 
ments such as the concentration of the cluster galaxy distribution or the dominance of the BCG as 
a possible second parameter to further refine the mass measurements. 

Figure [2] shows the color-magnitude distributions of galaxies with r/r 2 oo < 0.5 at z ~ 0.6. 
We have stacked 203 clusters with 0.55 < z p h < 0.65 to obtain the observed distribution (left 
panel). The red sequence is clearly visible. The patch-weighted background distribution is shown 
next in the middle panel. The distribution after background subtraction is shown in the right panel. 
A prominent red sequence appears after the subtraction. The solid line is the fitted red sequence at 
the effective redshift. The solid dots indicate the model M* at this redshift and the dashed vertical 
line indicates 1.5 magnitude fainter than M*. 



3.2. Semi-Empirical k + e corrections 

3.2.1. Red Sequence k + e models 



When computing statistics like red (or blue) fractions across a range of redshifts, the use of 
an evolving luminosity limit is crucial since the a bsolute luminosity of gala xies changes by about 
one magnitude from z = to z ~ 1. The original iButcher & Oemlen (|1984l) prescription of using 
a fixed rest-frame My = —20 does not produce a consistent galaxy population across a range 
of redshifts. This bias is especially severe for studies which span a large redshift range. The 
known correlation between luminosity and color can also mimic the evolution of the blue fraction. 
However, there is considerable unce rtain ty in the measu rement of luminosity evolution beyond 
z = 0.6. Work bv lFaberetallJ2005l) and bell etall d2004h showed that DEEP2 and COMBO-17 
surveys yielded very si milar results out to z ~ 1 with respect to the £ ?-band luminosity density for 
red galaxies. However, iBrown et al.1 (|2007|) and IScarlata et al.1 (120071) found a milder luminosity 
density evolution in the NOAO Deep Wide-Field and COSMOS surveys. 

The luminosity function was determined empirically from background- subtracted galaxy counts 
in the red seque nce of the clus t ers. A more detailed study of the RCS cluster luminosity function 
can be found in iGilbank et al.1 (|2008l) : our assumed luminosity function is within 0.1 mag of the 
values listed in their Table I. These values of M* (in the z' bandpass) as a function of redshift 
are consistent with the combination of /c-correction for an earl y-type galaxy spectrum along with 
passive evolution from a z ~ 3 single starburst population (see lGilbank et al.1 for details). 
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While our exact value of M* may differ slightly from that measured by others, the evolution 
in M* within our own sample is robust and will not influence trends with redshift. We explored the 
possible effects of uncertainty in the cluster galaxy luminosity function. A systematic change of 
±0.25 mag per redshift in M* for the cluster galaxies produced a change in red fraction less than 
la. 

Our goal is to obtain a consistent rest-frame flux for the cluster galaxies across the redshift 
range we are probing. We first use our model M* to convert the apparent z' magnitudes of red 
sequence galaxies to evolution-corrected luminosities at the redshift of the cluster. The existence 
of the red sequence then allows us to match the color-luminosity distribution across redshift. The 
observed slope can be used to change the apparent R — z' color i nto units of C* (R — z' color at 



M*), without appealing to model dependent color evolution (e.g.. iBell et al]|2004l : [Loh & Strauss 



20061). Figure [3] shows the color-luminosity distribution for two redshift slices: 0.45 < z p h < 0.55 
(first panel from the left) and 0.85 < z v ^ < 0.90 (3rd panel). The axes are in units of M* and C*, 
the color relative to the red sequence. 



3.2.2. Differential k-correction 

While model independent, the approach from the above section does not take into account 
the differential /^-corrections between galaxies of different spectral types for a fixed observed z' 
magnitude. To this end, we use the observed color blue-ward of the red-sequence to infer the 
relative /c-correction beyond the flux relative to a M * at that epoch , which gives a comparable 



k + e-correction to more traditional approaches (e.g..: lYee et al.ll2005|) 



The left panel of Figure |4] shows the synthetic R c — z' color for galaxies of diffe r ent sp ectral 



types computed using real galaxy template spectra from iColeman. Wu & Weedmanl (|1980l) aug- 
mented with two starburst models. The empirically determined red-sequence color as a function of 
redshift is shown in red. Since the red sequence includes passive evolution of the stellar population, 
it does not match the non-evolving colors derived from local ellipticals at high redshift. Plotted 
on the right panel of Figure |4] are the synthetic colors, relative to the local ellipticals, for galaxies 
of different spectral type as a function of redshift. The y-axis can be compared directly with the 
y-axis of the red-sequence corrected color-magnitude distribution of Figure [3] From the relative 
color of each galaxy, we assign a spectral type to each, interpolating if the color is in between 
types. The relative fc-correction in z'-band, used to shift a given galaxy to a lower intrinsic lumi- 
nosity, is inferred from the differential fc-correction relative to the fiducial red sequence template 
at each reshift. 



Figure [3] shows the changes in the color-magnitude distributions before (1st and 3rd panels) 
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and after (2nd and 4th panels) applying the relative ^-correction to the z'-band flux. This additional 
relative /c-correction corrects blue galaxies for on-going star-formation by pushing them to fainter 
magnitudes, such that the bluest faint galaxies, which would otherwise just make our luminosity 
cut, are pushed below the magnitude limit and are so removed from our sample. This effect is 
clearly visible at higher redshift by examining blue galaxies in the 3rd and 4th panels of Figure 
[3] (0.85 < z p h < 0.90). The numerous blue galaxies brighter than M* are corrected to fainter 
magnitudes but are still bright enough to be included in the sample. If we neglect the differential 
evolution between red sequence and blue galaxies, the derived red fraction at r 2 oo is reduced by 
2% at z ~ 0.5 and by 14% at z ~ 0.8. 



3.3. Cluster Red Fractions 



There are a number of different ways in which red or blue galaxy fractions are defined in the 
litera ture. Some use a fixed rest- frame color cut, say rest-frame B — V < 0.8 mag for blue galaxies 
(e.g. jRakos & Schombertll 19951) . Definitions like these usually yield a large number density evo- 
lution of early-t ype galaxies that are grossly inconsistent with the hiearchical model of structure 
formation (e.g.. IWolf et alj|2003|) . Others, such as the original iButcher & Oemlen (U.984) studies, 
use a fixed co lor width A fro m the observed re d sequen c e, say , A(B — V) = 0.2 mag for red 
galaxies (e.g.. bell et al.ll2004|) . As reviewed by lAndreonl (|2005r) . a constant width bracketing the 
red sequence is not optimal because galaxies of different spectral types drift in and out of the fixed 
region A for even the most generic models of stellar population. Claims of the evolution the blue 
(or red) fraction may just be a consequence of systematic spectral types drifting in and out of the 
width at different redshifts. 

Here, we use an empirical approach that attempts to separate the red sequence galaxies from 
the rest of the population, and divide the remaining galaxies amongst different fixed spectral types. 
We first construct the k + e corrected R c — z' color distribution by summing over the counts from 
all galaxies brighter than M* + 1.5. The distribution for redshifts slices of 0.45 < z p h < 0.55, 
0.65 < z p h < 0.75 and 0.85 < z ph < 0.90 are shown on the bottom panels of Figure [5J We 
then fit the red side (split by the red sequence peak) of the distribution with a double Gaussian 
of a single mean, including a narrow peak and a broader wing. Because the intrinsic spread 



of the red sequence is small (< 0.05 maa: iBower. Lucev & Ellisl Il992l : iBlakeslee et al.l 12003 
Lopez-Cruz. Barkhouse & YeeEo04l : Icool et al 2006 : lMei et al.ll2006|) and there should be on av 



erage no excess of galaxies with redder colors associated with these clusters, this half of the color 
profile should be characteristic of the color uncertainty in our observations and should be symmet- 
ric about the red sequence peak. This peak is prominent enough at all observed redshifts that red 
sequence galaxies are thus robustly isolated. The red fraction, f R ■ = 1 — f B , the complement of the 
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traditionally measured blue fraction, is estimated by mirroring the double Gaussian about the 
mean model distribution, and normalizing by the total galaxy counts (red and b lue, smoothed with 
a spli ne curve). We note here that our approach is in spirit similar to the original lButcher & Oemler 
(|1984|) studies where the color width of their re d sequence galax ies was indeed « 0.2 mag. This 
approach is also similar to the one suggested by Andreonl (|2005l) of using the valley between the 



two distributions as a divider (e.g. Jstrateva et al 



20011) and more r ecent bi-gaussian fitting to the 
local field color distribution (at fixed luminosity) by iBaldry et al. 1 d2004h . Our approach has the 
added advantage of estimating a more robust red fraction in the absence of a well-defined valley 
(e.g., the first distribution shown in Figure 0), and when there is a departure from bi-gaussianity. 
Furthermore, the observed red distribution seems to broaden at increasing redshift. Our approach 
automatically captures this broadening, presumably due to the larger photometric uncertainties 
of the individual sources at fainter apparent magnitudes. Both the valley di vision and the fixed 
width approach w ill not take this into account, creating a susceptibility to an lEddingtonl -like bias 
(Ed dingtonlll913l) which might mimic the redshift evolutionary trend we are trying to detect. 



3.4. Cluster-to-Cluster Variation 

It has been noted from many previous studies of the Butcher- Oemler effect that even at fixed 
redshift, there is considerable variation in blue fraction. Some of this variation is due to the het- 
erogeneous manner in which the clusters were selected, while others are intrinsic in the sense that 
even at a given epoch, clusters have various physical morphologies (e.g., Bautz-Morgan Class) 
and states of evolution. We try to capture this variation between clusters by performing a jackknife 
analysis on our stacked clusters at a given redshift slice. This is done by dividing the cluster sample 
further into 15 to 20 subsamples, depending on the number of clusters involved. Each subsample is 
then removed, and the same red fraction analysis is performed. The error-bars displayed in Figure 
|6] are obtained by this approach. We note here that this approach also captures other variations 
within our data set, like small residual differences in patch to patch relative photometry. 



4. Results and Discussion 

4.1. The Butcher-Oemler Effect 

The cumulative red fraction as a function of redshift is shown in Figure [6] The fraction of red 
galaxies decreases as a function of redshift for all scaled radii r/r2oo- Particularly relevant to the 
Butcher-Oemler phenomenon is the gradual decline in red fraction at the cluster core (0.25 r 2 oo) 
from ~ 95% at z ~ 0.5 to ~ 80% at z ~ 0.9. The red fraction value of our lowest redshift bin 
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is high compared with the original iButcher & Oemlerl (|l984|) study and many subsequent analyses 
which showed red fractions of about 0.7 at z ~ 0.4. Part of this discrepancy comes from the way in 
which the red fraction is measured; using a fixed color cut in the bimodal distribution will tend to 
increase the blue fraction if the observational scatter is significant. Our values are mo re consistent 
with recent studies using X-ray- selected cluster samples at z ~ 0.3 (|Wake et al.ll2005l) and z ~ 0.6 
(|Andreon. Lobo & Iovinoll2004l) . 



At face value, our result is inconsistent with the K-band analysis of iDe Propris et al.l (|2003|) 
which probed a similar redshift range, but was based on a more massive (on average) and het- 
erogeneous sample of 33 optical, X-ray and radio-selected clusters. They fitted their data with a 
constant /# ~ 0.1 out to z ~ 0.9. They confirmed the Butcher-Oemler effect as observed in an 
optically-defined sample, but noted little evolution when using a K-band selected sample. They 
concluded that any evolutionary effect may be due to a population of bursting dwarf galaxies, 
whose K-band luminosities (hence the stellar mass of the galaxies) are too faint for inclusion in a 
consistent K-band galaxy selection. 

Our use of k + e corrected magnitude in the ^'-bandpass is not a perfect proxy for galaxy 
stellar mass but is less affected by recent star formation and has a much smaller k + e correction 
compared with other optical bandpasses. Hence, we expect a reasonable agreement with K-band 
results. Howe ver, we find that ou r analy s es differ significa ntly in the division of red versus blue 
galaxies. Like butcher & Oemlerl (|1984|) . IDe Propris et al.1 use a fixed A(B — V) — 0.2 mag but 
fc-corrected to their observed bandpasses using an eff ective SEP derived from Sa + Sc mixture to 
give an effective A(color) at higher redshift. Unlike butcher & Oemlerl whose color histograms 
have an observed red sequence width of ~ 0.2 mag, the color histograms from IPe Propris et al. 



display a narrower red sequence width compared with the effective A(color) - often as large as 
0.5 mag - used in their analysis. Their chosen A(color) allows the inclusion of galaxies of later 
spectral type in the red fraction, suppressing their blue fraction measurements relative to ours, 
which model the observational scatter. On examination of the individual color histograms of their 
high redshift clusters, we find that our method would produce substantial blue counts beyond the 
red sequence, and hence we find no evidence of inconsistency between our results and theirs. 



We further note that the IPe Propris et all s results are derived using a fixed aperture of rela- 
tively small size (0.25 — 0.45/i" 1 Mpc), whereas our red fractions are measured using apertures of 
scaled radii r/r 2 oo- Observations made within a fixed metric radius tend to sample only the redder 
inner cores for massive clusters, and at z > 0.6, their sample is likely biased toward the most 
massive systems for which the r 2 oo values are typically as large as 1.5 to 2.0h~ 1 Mpc. As Figure 6 
for the RCS sample shown, at radii smaller than 0.5r 2 oo> there is only a moderate evolution in the 
red fraction between z ~ 0.4 and 0.9. 



Our definition of red fraction is sensitive to our assumption that the red sequence color spread 
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is symmetric, and dominated by observational uncertainties. This appears to be the case for most 
of our data. However, in our two lowest redshift slices 0.45 < z ph < 0.50 and 0.50 < z ph < 0.55, 
we find that the chosen filters are not ideal separators of the red and blue distributions. In these two 
bins, the blue galaxies scatter underneath the red sequence so that our reflection of the red wing 
of the distribution will tend to oversubtract by a small amount. This effect is more problematic at 
fainter magnitudes because of the larger fraction of blue galaxies. For this redshift bin we fit the red 
sequence separately in four different magnitude bins to better separate the two distributions. Even 
so, the red fraction measured in this bin may be slightly too large by at most 5%. The problem 
only occurs in the lowest redshift slice, since at higher redshifts, the blue cloud is more widely 
separated in color space from the red peak for these filters. Inclusion of data from bluer filters will 
allow a bette r separation of th e red and blue galaxies at the lower redshifts. An analysis using four 
color data (|Hsieh et al.ll2005|) over a smaller area will be explored in a future paper. 



4.2. Radial Dependence 



A population gradient is seen for all redshifts in that the red fraction is larger in the cluster 
cores at all redshifts. In Figure [6] we plot the red fraction as a function of redshift for galaxies 
within r/r 200 < 0.25, 0.25 < r/r 200 < 0.50, 0.50 < r/r 200 < 1.00 and 1.00 < r/r - mn < 2.00 



This result is similar to man y low-redshift stu dies of clus t ers (e .g. jAndreon et al.H2006l : lGoto et al 
20041 : lEUingson et alJ200lh . In their analysis,|GotoetaD (|2004|) found a break at r /r 200 ~ 0.3-0.8 



at their highest redshift bin (0.2 < z ph < 0.3). Our results indicate a similar change in the radial 
gradient of the red fraction at r/r 20 o > 0.50. We caution such interpretations in light of the 
somewhat uncertain mass estimates of the clusters. Both statistical and systematic uncertainties 
in the cluster masses may affect our measures of the evolution of the red fraction and its radial 
dependence. Mass estimates from B gc measurements have ~ 60% uncertainties for individual 
systems, which translate to 30% uncertainties in r 20 o- Population gradients will be smoothed on 
this scale, possibly flattening any radial gradients on this scale, but most likely not affecting the 
comparison of evolution between 0.25 r 2 oo and 1.00 r 2 oo- 

The decline in red fraction with redshift is significantly steeper in the outskirts (1.00 < 
^/^2oo < 2. 00) compared to the c ore ( r/r 2 nn < 0.25). Our resu lts qualitatively agree with ear- 
lier work by Ellingson et al.1 (|200l|) and lKodama & Bowed (|200ll) who found a milder decrease in 
red fraction at the core compared to the outskirts/field using an X-ray-selected clusters sample with 
spectroscopic redshifts. A possible interpretation is that this effect comes from a relative decline 
in the infall rates of bluer galaxies into clusters, as new additions to the cluster should generally 
lie at larger radii. We note that while systematic errors in our estimates of cluster masses, galaxy 
luminosity functions and /^-corrections may affect our calculation of red fractions, this differential 
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effect with respect to radius will remain, suggesting that the population in the cores of clusters 
changes with time differently from the population in the outskirts. 

Our result on the radial dependence of the red fraction requires an accurate determination 
of each cluster center. A detailed description for determining RCS cluster centers was given in 
Gladders & Yed (|2005|) . In brief, the location of a cluster is first selected from the color- position 
data-cube with a magnitude weighting towards galaxies brighter than M*; the positioning is then 
refined by an interative Gaussian kernel of 250/iT 1 kpc. To the extent that the RCS algorithm is 
accurate, then the centering should be better than < 250/2T 1 kpc. Note that clusters in reality are 
not round objects and often there are double clusters, or galaxy overdensities along filamentary 
structures. If a cluster is indeed elliptical and has a well-defined center, we argue here that we 
will get the right center to better than half our searching scale (i.e. 125/i -1 kpc). If the cluster 
or galaxy overdensity is filamentary, then our assumption of self- similar spherical clusters breaks 
down. Indeed our results at r < 0.25r 2 oo maybe affected by centering problems, but not for the 
larger radii. If we are not systematically deviating in our cluster center determination, then the 
jackknife procedure in our error calculation would capture the diversity and uncertainty of our red 
fraction determination. 

Because cluster galaxy populations show such marked radial gradients, any systematic varia- 
tions of the B gc - r 2 oo relation with redshift may also change our red fraction as a function redshift, 
with any systematic tendency to overestimate cluster masses leading to a trend towards a lower red 
fraction. However, a significant systematic is required: even if our estimates of r 20 o are off by a fac- 
tor of 3 at z ~ 0.8, compared to z ~ 0.5, we would still see an evolution with redshift. Currently, 
no such redshift dependence is found in the RCS mass estimates when the same cluster s ample is 
analysed using a global self-calibration pro cedure (e.g., Hull2003 : Majumdar & Mohrll2004 ) which 
reflects standard cosmological parameters (IGladders et al.ll2007l) . Ongoing wor k is currently bein g 
done to calibrate the mass of a subset of high redshift clusters using X-ray (IHicks et al.ll2007|) . 
multi-object spectroscopy (|Gilbank et al.l 120071 : Barrientos et al., in prep), weak lensing analyses 
(Hoekstra et al., in prep.) and mid-IR (Ellingson et al., in prep.) observations to constrain any such 
biases. 



4.3. The Blue Galaxy Distribution 

Figure [7] shows the blue distribution of galaxies after removing the fitted red distribution. The 
shaded vertical bands are the expected color for galaxies of different (late) spectral type, ranging 
from Sab (far right) to starburst (SB; far left), synthesized using a single representative SED for 
each type (cf. Figure H]). The width of the band reflects the finite redshift range in each subsample 
and the color variation as the filters shift along the SED with redshift. The blue distribution, 
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when compared with spectral type, appears to change slightly with redshift, becoming bluer at 
higher redshifts. The broadening towards the red tail at the highest redshift sample may be due 
to photometric uncertainty (0.10 — 0.15 mag), but this cannot account for the longer tail to the 
blue. We find that the median galaxy spectral type for non-red- sequence galaxies changes by 
approximately half a spectral type between z = 0.4 and z = 0.9, (dashed line in FigurefTl and the 
series of six blue squares in the middle panel of Figure lU, and the fraction of blue galaxies with 
colors bluer than Im increases from 10% to 40 %, between redshifts 0.5 and 0. 9. Studies of color 



bimodality in field galaxies (IStrateva et al.ll200ll : iBell et al.ll2004l : lBlantonll2Q06h suggest that field 
galaxies have colors consistent with slightly later types at these redshifts. 

The Butcher- Oemler effect has been attributed by some studies to a population of very blue 



starb urst galaxies in clusters at high redshift (e.g., iRakos & Schombertl 1 19951 : iDe Propris et al 



20031). We investigate this possibilty by examining the redshift behavior of the fractions of blue 
galaxies with colors bluer then Im spectral type. The contribution from galaxies with spectral 
type similar to unobscured starbursts (we take the color mid-point between the the two models) 
is < 5%, even for the highest redshift slice. This is an upper bound because we expect leakages 
from the more abundant redder galaxies to the blue tail due to photometric errors. Looking fainter 
into the luminosity function (> M* + 1.5) would increase the contribution from galaxies with later 
spectral types, but we show that the Butcher-Oemler effect as detected here is not driven by low- 
mass starbursts, and includes a component of brighter galaxies w ith c olors consistent with n ormal 
coeval field galaxies. Recently, works from Elbaz et al.1 (|2007r) and ICooper et al.1 (|2007br) have 
found an increase in the number of bright blue, normal (i.e. not low mass or starburst) galaxies 
in high density environments at z = 1 relative to z = (which they suggest may drive an inver- 
sion of the star-formation - density relation). Although our survey only probes rich clusters, our 
results suggest that the Butcher-Oemler effect is driven by normal galaxies, and are qualitatively 
consistent. 



4.4. Selection Effects from the Red Sequence Method 



By construction, our clusters are selected based on the red sequence method using z' and 
R c bandpasses, which will create a bias against finding clusters with very blue populations. Our 
measurement of the Butcher-Oemler effect may thus be considered to be a lower limit. However, 
we do not expect that such incompleteness will cre ate a large addi tional effect. Simulations of 
cluster detection efficiency using the RCS method (Gladdersl 120021) show that the least massive 
clusters in our sample (corresponding to B gc ps 300 (h§Q Mpc) 177 ) at z ~ 0.9 are detected at high 
efficiency with red fractions as low as 0.55. If the red fraction dips to only 0.2, this efficiency drops 
to 50%. More massive clusters may be detected at 80% efficiency at a red fraction of 0.2. Since we 
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do not calculate red fractions for individual clusters here, we cannot use this to create a quantitative 
correction. However, these efficiencies are still significant even in the worst case scenarios, and 
average values lie significantly far away from the red fractions where we expect a selection bias to 
occur. Moreover, no other method for finding high redshift clusters has identified a sample which 
contains a large fraction of massive clusters with these very blue colors, and the overall number 
of clusters found in the RCS survey appears to be consistent with current cosmologies. Thus, we 
conclude that the problem of missing very blue clusters of galaxies is not overwhelming, and the 
Butcher-Oemler effect is not likely to be significantly stronger than observed. Because of our red- 
sequence selection criterion, there is a possibility that our sample may have a slight relative bias 
toward the redder systems at higher redshift. If this is the case, then the observed Butcher-Oemler 
effect would be a slight underestimate of the true evolution. 



5. Summary and Conclusions 



We use Rc and z' observations of a sample of ~ 1000 clusters with mas ses greater than 
10 14 /i _1 M and 0.45 < z ph < 0.90 from the Red-sequence Cluster Survey (|Gladders & Yee 



2005|) to investigate the evolution of the color-magnitude relationship for cluster galaxies. We con- 
struct composite color-magnitude diagrams for a a number of redshift bins and ranges of cluster- 
centric radii, using cluster photometric redshifts and a statistical background subtraction technique. 
We correct each galaxy color distribution for an empirically-determined estimate of the cluster lu- 
minosity function and /^-corrections as a function of galaxy color, and investigate the bimodal 
distribution of excess galaxies in cluster fields. The large dataset and composite technique allows 
us to avoid many of the difficulties of previous investigations, where the luminosity and color evo- 
lution of different galaxy types and observational scatter blur the definitions of "red" and "blue" 
galaxies. Our analy sis is informed by recent s tudies that the red sequence of galaxies is intrinsi- 
cally narrow (e.g. jBower. Lucey & EUislll992|) . allowing us to fit an empirical description of the 
effects of observational scatter in our data, and separate galaxies in the red sequence from the rest 
of the cluster population. 

We find a decrease in the fraction of red-sequence galaxies (increase in the blue fraction) 
in clusters across this redshift range, qualitatively consistent with the Butcher-Oemler effect. Our 
measurements are difficult to compare quantitatively with previous measures of the blue fraction, 
and tend to yield somewhat higher red fractions because of our correction for scatter in the red 
sequence and explicit k-corrections for b lue galaxies. Because of this, our results indicate a m ilder 
evolution than some earlier studies (e.g., Butcher & Oemlei 1984 : lRakos & Schombert 1995 ). but 
a stronger evolution than that suggested by iDe Propris et alj (|2003|) using a small hetergeneous 
sample with K-band photometry. Our results are more consistent with recent studies of X-ray 
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selected cluster samples (e.g.. IWake et al]|2005|) . We argue that while the red sequence technique 
may miss some fraction of very poor, blue clusters, this effect is minor, and the Butcher-Oemler 
effect is unlikely to be significantly stronger than seen here. We see consistent radial gradients in 
clusters, in that red fractions are higher in clusters cores than within larger radii at all redshifts 
studied. We also find a mild differential in the rate at which populations evolve - more slowly 
within the inner cores of the clusters (r < 0.25 r 2 oo) tha n in the outer regions (1.00 < r/r 2 nn < 
2.00) , consistent with previ ous results at 0.3 < z < 0.6 (lEUingson et al.ll200ll : iKodama & Bower 
200 ll : lAndreon et al.l |2006|) . This differential may indicate that cosmologically-driven infall of 
bluer galaxies into clusters may be partially responsible for the evolution in populations, as these 
galaxies will preferentially be found at larger radii. The decline in i nfall rates in low -density 
cosmologies then contributes to the decline in this population over time (|Ellingsonll2003|) . 



The color distribution of non-red- sequence galaxies appears to be consistent with a popula- 
tion of normal spirals, and shows a gradual evolution in the median color with redshift. There 
is only a minor contribution from galaxies with colors similar to unobscured starbursts. This 
places robust limits on interpretations of the Butcher-Oemler effect as being driven by starburst 



ing dwarf galaxies that subsequently fade dramatically in clusters (e.g jRakos & Schombertlll995 



De Propris et al.l 120031) . However, the two-color data used here is insufficient to address quan- 
titatively the star formation rates in cluster galaxies, or assess the possibility that obscured star- 
burst and post- starbur st galaxies play an important role in the evolution of cluster populations 
(|Poggianti et al.lll999l) . Additional studies of galaxy spectral energy distributions, morphologies 
and star formation rates are underway to trace the evolution of individual galaxies in clusters in 
this sample. 
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Table 1: 



red shift range 


Number of Clusters 


0.45 < z < 0.55 


171 


0.55 < z < 0.65 


203 


0.65 < z < 0.75 


229 


0.75 < z < 0.85 


263 


0.85 < z < 0.90 


118 




0.5 0.6 0.7 0.8 0.9 

Redshift 



Fig. 1. — The redshift histogram of RCS clusters used in our analysis. The three segments indicate 
richness bins of B gc > 800 (grey for rich clusters), 500 < B gc < 800 (hatched) and 300 < B gc < 
500 (white for poor clusters). 
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Fig. 2. — (Left) The observed color-magnitude distribution for all non-stellar sources with 
r < 0.5 r 20 o for 203 clusters with 0.55 < z p h < 0.65. (Middle) Patch- weighted background 
distribution (see text). (Right) Background-subtracted distribution of galaxies, calculated from the 
two preceding panels. The white dot indicates M* at the median redshift (~ 0.61) while the dashed 
vertical line indicates M* + 1.5, the consistent flux limit used for the red fraction analysis. The 
completeness of the photometry for extended sources (z' < 23.1 and R c < 24.2) lies beyond the 
boundary of these figures. 





Fig. 3. — Color-magnitude distribution after k + e correction for 0.45 < z ph < 0.55 (1st and 2nd 
panel from the left) and 0.85 < z p ^ < 0.90 (3rd and 4th). The axes of the plots are in units of 
M* and C* at their respective redshifts. In the 1st and 3rd panel, the z'-band flux is corrected 
consistently for red sequence galaxies only. In the 2nd and 4th panel, a differential ^-correction is 
based on the expected color of galaxies of different spectral type. This latter correction not only 
tightens the red sequence at low redshift, but also substantially reduces the counts of blue M < M* 
galaxies at high redshift. 
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Fig. 4. — (Left) Synthetic R r — z' colors for galaxies of different spectral type computed using 
galaxy template spectra from Coleman. Wu & Weedman (|l98ol) and model starburst (SB1, SB2) 
galaxy spectra. The evolving red-sequence M* model is shown in red. The series of squares are the 
median color of the non-red- sequence galaxies from RCS clusters (see section H31 and Figure|7]) as 
a function of cluster redshift. (Right) Synthetic colors relative to a fiducial ~ M* elliptical galaxy 
as a function of redshift. The y-axis is in units of C*, the color relative to the colors of ellipticals 
from the CWW models. The vertical dashed lines bracket the redshift range of our sample. 
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R - z (observed) R - z (observed) R - z (observed) 

Fig. 5. — The top panels are the background subtracted color-magnitude distributions. The bottom 
panels are the color distribution after k + e correction. The histogram and the solid line shows the 
distribution after applying a differential ^-correction for non-red sequence galaxies (see section 
13.2.21) while the dashed lines show the distribution before such a correction was made. The red 
curve is our model for the red distribution used to compute the red fraction. The plots for are 
redshift slices: 0.45 < z ph < 0.55 (left), 0.65 < z ph < 0.75 (middle) and 0.85 < z ph < 0.90 
(right). 



-26- 



1.0 




0.4 0.5 0.6 0.7 0.8 0.9 

Redshift 



Fig. 6. — The red fraction as a function of redshift for galaxies with r/r 2 oo < 0.25 (core), 0.25 < 
r/r 2 oo < 0.5, 0.5 < r/r 2 oo < 1-0 and 1.0 < r/r 2 oo < 2.0 (outskirts) and brighter than M* + 1.5. 
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Fig. 7. — The blue distributions of galaxies after correction for the red sequence (solid histograms). 
The shaded bands in each plot are the expected color for different spectral types, ranging from Sab 
(far right), Sbc, Im, and starburst (SB1, SB2; far left), synthesized using a single representative 
SED for each type (see text). The width of the band reflects the finite redshift range in each 
subsample and the color variation as the filters shifts along the SED with redshift. The plots are 
for redshift slices: 0.45 < z ph < 0.55 (left), 0.65 < z ph < 0.75 (middle) and 0.85 < z ph < 0.90 
(right). The vertical dashed lines indicate the median colors of the distributions, also plotted in 
Figure 4. 



